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Solids with spatial variations in the crystalline axes naturally evolve into cells or grains separated 
by sharp walls. Such variations are mathematically described using the Nye dislocation density 
tensor. At high temperatures, polycrystalline grains form from the melt and coarsen with time: the 
dislocations can both climb and glide. At low temperatures under shear the dislocations (which 
allow only glide) form into cell structures. While both the microscopic laws of dislocation motion 
and the macroscopic laws of coarsening and plastic deformation are well studied, we hitherto have 
had no simple, continuum explanation for the evolution of dislocations into sharp walls. We present 
here a mesoscale theory of dislocation motion. It provides a quantitative description of deformation 
and rotation, grounded in a microscopic order parameter field exhibiting the topologically conserved 
quantities. The topological current of the Nye dislocation density tensor is derived from a micro- 
scopic theory of glide driven by Peach-Koehler forces between dislocations using a simple closure 
approximation. The resulting theory is shown to form sharp dislocation walls in finite time, both 
with and without dislocation climb. 

PACS numbers: 61.72.Bb,61.72.Lk,62.20Fc 



Crystals, when formed or deformed, relax by develop- 
ing walls. Common metals (coins, silverware) are poly- 
crystalline: the atoms locally arrange into grains each 
with a specific crystalline lattice orientation, separated 
by sharp, flat walls called grain boundaries. When met- 
als are deformed (pounded or permanently bent) new 
cell walls form inside each grairi^. Until now, our only 
convincing understanding of why crystals form walls has 
been detailed and microscopic. Our new theory provides 
an elegant, continuum description of cell wall formation 
as the development of a shock front - a phenomenon hith- 
erto associated with traffic jams and sonic booms. 

Continuum theories of dislocation dynamics are not 
neTO 4 i 5 4 6 4 7 Various approaches have been used to study 
dislocation pattern formation. Simplified diffusion mod- 
els have been used to describe persistent slip bands&, dis- 
location cell structures during multiple slipS, and dislo- 
cation vein structures^. Models similar in spirit to ours 
have recently been analyzed both for single-slipii*i2ii£ 
and in three dimensions^. None of these approaches 
have yielded the sharp, wall-like singularities character- 
istic of grain boundaries and deformation-induced cell 
structures (see, however—). 

Continuum plasticity theory, as used in practical engi- 
neering codes, describes the deformation gradient (diUj, 
for displacement field u) as the sum diUj = (3fj +/? P - of an 
elastic, reversible distortion /3 E and a permanent plastic 
distortion /3 P . The plastic distortion involves the cre- 
ation of microscopic dislocations, and cannot be written 
as a gradient of a single-valued displacement field. Inte- 
grating around a loop L enclosing a surface S, the change 
in such a hypothetical plastic distortion field Au p can be 
written using Stokes theorem as 
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where the Nye dislocation density tensor 
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measures the net flux of dislocations a, tangent to t, with 
Burgers vector b, in the (coarse-grained) neighborhood 
of x. The microscopic statement that dislocations cannot 
end implies dtpij = 0, so the time evolution must be given 



in terms of a currentSii& J: dpij/dt ■ 
From equation J2J we see thatii 



Ja = opudt. 
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Thus the natural physicist's order parameter (the topo- 
logically conserved dislocation density p) is a curl of the 
common engineering state variable (the plastic distortion 
field /3 P ). The focus of our manuscript will be the deriva- 
tion of an evolution law Q appropriate for scales large 
compared to the atoms but small compared to the cell 
structures and grain boundaries. 

Notice that we only consider the net density of disloca- 
tions. We ignore the geometrically unneccessary disloca- 
tions (those with opposing Burgers vectors which cancel 
out in the net dislocation density) because they do not af- 
fect the long range strain fields or the misorientations at 
grain boundaries and cell walls. Macroscopically they are 
known to dominate dislocation entanglement and work 
hardening, and are included in previous continuum the- 
ories of plasticitj&iiiiSiiiiSiiS. Much of the macroscopic 
cancellation in net dislocation density comes from the 
near alternation of the net rotations in the series of cell 
walls^S. Our focus on the sub-cellular, subgrain length 
scales and our current omission of dislocation tangling 
make keeping only the net dislocation density natural 
for our purposes. We also do not explicitly incorporate a 
yield surface, because we hope eventually to explain work 
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hardening and yield surfaces as properties which emerge 
from the intermediate length-scale theory. 

We motivate our evolution law from the microscopic 
Peach-Kochler force on a section of a dislocation line 
due to the stress field a present at that point: f PK = 
—EijktjbiCTki- The current J of a single dislocation mov- 
ing with velocity v is Jy = £u m tibjV m 5(£). If we treat 
dislocation glide and climb on an equal footing (a crude 
model for metals at high temperatures), the dislocation 
will move in the direction of the applied force v oc f PK . 
For a single dislocation, tuning the climb component with 
a parameter A, we find 
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Here D is a materials constant with units of 
[length] 2 [time] /[mass] giving the mobility of dislocation 
glide. At A = climb and glide have equal mobilities, and 
at A = 1 J is traceless and hence only glide is allowed. 
Coarse graining over many dislocations, 
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)] is an eighth rank 
isotropic tensor specific to Peach-Koehler model, and 
Pkmrs i s a hig ncr order Nye dislocation tensor. Thus that 
the evolution of p depends upon the new object p( 4 \ 

We now need to perform a closure approximation, writ- 
ing p( 4 ) in terms of p (as in Hartree-Fock and in theories 
of turbulence). The only choice p^ — > p® p that guar- 
antees a decrease of elastic energy with time^l is 

^tf&iWCn (6) 
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where C has units of distance. The resulting evolution 
law for the plastic distortion tensor is 
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where p is a curl of f3 p (equation E}, and the stress <r 
(due to the long-range fields of the other dislocations) is 
written in Fourier space as^ 
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For simplicity, we present solutions to our three- 
dimensional theory for the special case of systems where 
the dislocation density varies only along one and two di- 
mensions. In one dimension with variations only along 
z, our theory becomes local, with stress 



tr(z) = -p 
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Nye dislocation tensor 

Pxj=d Z Pyj, P yj = -d z Pl 



and evolution law 
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with 



F((3 P ) = 
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+ ^ P J+(3 P y 2 )+j^((3 P x + (3 P y j 2 }. (12) 

Interestingly, w((3 p ) is the linear clastic energy density of 
the system, the effective force F((3 P ) is proportional to 
the ^-component of the Peach-Kochler force, and the rate 
of change of the elastic energy density is simply —F 2 ((3 P ) 
(so the energy decreasing condition is explicit). 

Our equation for /3 P in one dimension is hyperbolic: 
the method of characteristics may be applied, where the 
characteristics are the parameterized curves (t(s),z(s)) 
with dt/ols = l,dz/ds = F{0 P ). Along these curves, 
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implying (from equation lllfl that /3 P is constant along 
the characteristics. As with many hyperbolic systems, 
we observe shock formuation in finite time. Our ana- 
lytic analysis of the shock asymptotics (in preparation) 
exhibits the jump in (3 P , wall in p, and stress jump de- 
scribed numerically below, with a l/yi time-dependence 
for the surrounding continuum (3 P . 

We simulate systems of spatial extent L in one or two 
dimensions, with periodic boundary conditions and no 
external stress. The initial plastic distortion field j3 p is 
a Gaussian random field with mean square amplitude /3 P 
and root-mean-square decay length approximately L/16. 
In one dimension, we use the upwind scheme 2 ^ as imple- 
mented by Press et aU^. In two dimensions, we calculate 
er using Fourier methods and regularize our equations by 
adding a fourth-order numerical viscosity — i4TV 4 /3 p - to 
equation [7| (A second-order viscosity V 2 /?? was not as 
successful in suppressing the instabilities.) In one dimen- 
sion, we have checked that the upwind scheme produces 
the same solution as the viscosity regularization. Unlike 
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shock waves in fluids, where conservation laws tell us that 
the upwind/ viscosity solution is the correct one, here the 
dynamics of the walls is not determined by the contin- 
uum theory. (We speculate that the wall dynamics reg- 
ularization may be non-trivial and stochastic, emerging 
out of avalanches and critical depinning phenomena on 
sub-cell length scales - observed both experimentally 25 
and theoretically 2 ^ in systems with only one active slip 
system.) 
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FIG. 1: Plastic distortion tensor component 0fj in one di- 
mension allowing only glide motion, after time t = 20L 2 /-D/x, 
with 2048 mesh points. The shocks or jumps in the values 
correspond to the cell walls. A stress-free wall (satisfying the 
Frank conditions) here would have no jump in 0^ x , 0y y , or 
P + F 

Figure 2] shows the final state of the plastic distor- 
tion /3 P in a one-dimensional system evolving under dis- 
location glide only (corresponding to plastic deformation, 
which experimentally leads to cell structure formation). 
Jump singularities in /3 P form after a short time, repre- 
senting a wall of dislocations and an abrupt change in lat- 
tice orientation. The cell walls in our model are not grain 
boundaries, because they do not correspond to pure, sim- 
ple rotations (they do not satisfy the Frank conditions 2 ^). 
There is a jump in the stress across each cell wall: the 
glide component of the forces from the two neighboring 
cells on p in a wall is equal, opposite, and compressive. 
Our model therefore predicts that the initial formation 
of cell walls during low-temperature plastic deformation 
is driven by a stress jump due to non-cancellation of the 
stress field of the constituent dislocations. As our system 
evolves, new cell walls form (intriguingly similar to the 
cell refinement seen experimentally^* 2 *.) and the existing 
cell walls evolve to reduce their stress jumps. 

Figure [3 shows the evolution of the plastic distortion 
in a system allowing both glide and climb, corresponding 
to high-temperature annealing of a plastically deformed 
structure. Here the initial singularities again arc not 
grain boundaries, but walls of dislocations with net stress 
jumps. These walls then move and coalesce in response 
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FIG. 2: The ^-component of the plastic distortion tensor 
in one dimension up to time t = 22L 2 /Dfi, with 2048 mesh 
points. The evolution allows both glide and climb motions. 
The walls move and coalesce until only a single wall survives. 



to the net force due to other dislocations, coarsening the 
resulting grain boundary structure. 




FIG. 3: The yz-component of the plastic distortion ten- 
sor allowing only glide, in two dimensions after time t = 
9.15L 2 /£>M- 

Figure [3] shows the final state of /3 P in a two- 
dimensional plasticity simulation (with glide and not 
climb). Notice the formation of sharp walls (mathe- 
matically demonstrated above only in one dimension). 
These walls separate relatively unstrained regions which 
we identify with the cells formed during plastic defor- 
mation. As in the one-dimensional simulations, the cell 
walls formed at short times have stress jumps. 

Generically we expect our model to form sharp walls 
in finite time. Special initial conditions, however, may 
not form walls. Indeed, the case of a single slip system 
does not form singularities within our model — perhaps 
explaining why cell wall formation is not observed in 
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Stage I hardening, or in dislocation-dynamics studies of 
single slip in two dimensions otherwise similar to our ap- 
proac h 1 ^ 12 ' 1 ? . 

Our dislocation dynamics theory is minimalist: it ig- 
nores many features (geometrically unnecessary disloca- 
tions^, slip systems^, dislocation tangling, yield sur- 
faces, nucleation of new dislocations) that are known to 
be macroscopically important in real materials. It does 
incorporate cleanly and microscopically the topological 
constraints, long-range forces and energetics driving the 
dislocation dynamics. As hypothesized by the LEDS 
(low-energy dislocation structures) approaches, we find 
that a dynamics driven by minimizing energy (omit- 
ting tangling and nucleation) still produces cell boundary 
structures. Finally our theory, to our surprise, initially 
forms sharp walls that are not the usual zero-stress grain 
boundaries. 

In condensed-matter physics, crystals are anomalous. 



Most phases (liquid crystals, supcrfluids, superconduc- 
tors, magnets) respond smoothly in space when strained. 
Crystals, when rotationally strained, form sharp walls 
separating cells or grains. Our analysis suggests that the 
formation of these walls can be understood as a shock 
formation in the dislocation dynamics. 
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